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Abstract 

This  paper  presents  the  numerical  modeling  and  optimization  of  a  spectrum  splitting  photovoltaic-thermoelectric  (PV-TE)  hybrid 
system.  In  this  work,  a  simulation  model  is  established  in  consideration  of  solar  concentration  levels  and  several  heat  dissipation  rates. 
Exemplarily,  the  performance  of  a  hybrid  system  composed  of  a  GaAs  solar  cell  and  a  skutterudites  CoSb3  solar  thermoelectric  generator 
(TEG)  is  simulated.  Analysis  under  different  conditions  has  been  carried  out  to  evaluate  the  electrical  and  thermal  performance  of  the 
hybrid  system.  Results  show  that  the  cutoff-wavelength  of  the  GaAs-CoSb3  hybrid  system  is  mainly  determined  by  the  band  gap  of  solar 
cell,  when  the  solar  concentration  ratio  is  ranged  between  550  to  770  and  heat  transfer  coefficient  h  =  3000-4500  W/m2  K,  the  hybrid 
system  has  good  electrical  performance  and  low  operating  temperatures.  Based  on  the  analysis  of  the  GaAs-CoSb3  hybrid  system,  guide¬ 
lines  for  the  PV-TE  system  design  are  proposed.  It  is  also  compared  with  a  PV-only  system  working  under  the  same  cooling  condition; 
results  show  that  the  PV-TE  hybrid  system  is  more  suitable  for  working  under  high  concentrations. 

©2012  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Hybrid  solar  energy  conversion  technologies  always 
focus  on  the  cascade  use  of  incident  broadband  solar  spec¬ 
trum.  By  combining  the  electrical,  thermal  and  chemical 
solar  energy  conversion  methods,  better  energy  conversion 
efficiency  and  output  power  can  be  achieved.  The  photovol¬ 
taic-thermoelectric  (PV-TE)  hybrid  system  is  one  of  the 
demonstrated  technologies  that  are  useful  for  harvesting 
solar  energy.  Current  research  directions  of  PV-TE  hybrid 
technology  are  focused  on  the  direct  production  of  electric¬ 
ity,  the  hybird  electric  vehicle  energy  saving  and  the  cooling 
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system  application.  Yu  et  al.  (2008),  Leonov  et  al.  (2010), 
Van  Sark  (2011),  Yang  and  Yin  (2011)  and  Wang  et  al. 
(2011)  attached  thermoelectric  generators  to  the  back  side 
of  PV  modules  for  energy  harvesting,  the  performance  of 
these  hybrid  systems  is  investigated  theoretically  and  exper¬ 
imentally  in  their  studies.  Zhang  et  al.  (2009),  Zhang  and 
Chau  (201  la, b)  and  Fan  et  al.  (2010)  implemented  PV-TE 
hybrid  systems  for  hybrid  electric  vehicles,  in  which  thermo¬ 
electric  generators  were  used  to  collect  waste  heat  of  internal 
combustion  engines  and  solar  cells  for  converting  solar 
energy,  their  researches  were  focused  on  power  control  strat¬ 
egies  of  the  PV-TE  hybrid  system.  And  in  the  study  of  Jradi 
et  al.  (2011),  an  integrated  solar-driven  thermoelectric 
cooling  system  was  developed  to  produce  fresh  water  from 
humid  air. 
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Nomenclature 

A 

area  (m2) 

T 

temperature  (K) 

A’1 

n-type  leg  cross-section  area  (m2) 

V 

y  oc 

open  circuit  voltage  of  solar  cell  (V) 

Ap 

p-type  leg  cross-section  area  (m2) 

y TE 

Seebeck  voltage  of  TE  module  (V) 

c 

speed  of  light  in  a  vacuum  (m/s) 

Cs 

geometrical  concentration  ratio 

Greek  letters 

CXmax  concentration  ratio  when  the  system  achieves 

oc 

Seebeck  coefficient  of  TE  material  (V/K) 

its  maximum  efficiency 

&HC 

absorptivity  of  heat  collector 

Eg 

solar  cell  band  gap  (eV) 

ttpy 

absorptivity  of  solar  cell 

Ein 

incident  solar  radiation  (W/m2) 

Pr, 

temperature  coefficient  of  efficiency  (K_1) 

EPV 

solar  energy  directed  to  CPV  subsystem  (W) 

Pise 

temperature  coefficient  of  short  circuit  current 

ETe 

solar  energy  directed  to  TEG  subsystem  (W) 

(K-1) 

EQE 

solar  cell  external  quantum  efficiency  (%) 

Pvoc 

temperature  coefficient  of  solar  cell  open  circuit 

FF 

solar  cell  fill  factor  (%) 

voltage  (K-1) 

m 

AM1.5D  spectral  irradiance  (W/m2  nm) 

y 

feasibility  of  system  working  under  concentra- 

hcool 

heat  transfer  coefficient  of  cooling  system  (W/ 

tion  (m2  K/W) 

m2  K) 

8 

emissivity 

h 

unc 

natural  convection  heat  transfer  coefficient  (W/ 

rj 

efficiency  (%) 

m2  K) 

Oopt 

efficiency  of  optical  system  (%) 

hP 

Planck  constant  (m2  kg/s) 

K 

thermal  conductivity  of  TE  material  (W/m  K) 

Isc 

short  circuit  current  of  solar  cell  (A) 

2 

wavelength  of  photons  (nm) 

We 

current  of  TE  module  (A) 

K 

cutoff  wavelength  (nm) 

J sc 

short  circuit  current  density  (A/m2) 

a 

electrical  conductivity  of  TE  material  (S/m) 

K 

thermal  conductivity  (W/m-K) 

&  SB 

Stefan-Boltzmann  constant  (W/m2  K4) 

kp 

Boltzmann  constant  (m2  kg/s2  K) 

KHc 

thermal  conductivity  of  heat  collector  (W/K) 

Superscripts 

L 

length  of  TE  module  legs  (m) 

P 

p-type-leg  of  TE  module 

n 

quantity  of  finite  elements  of  TE  module 

n 

n-type-leg  of  TE  module 

nf 

diode  ideality  factor  of  solar  cell 

P 

generated  power  (W) 

Subscripts 

q 

elementary  charge  (C) 

AIR 

ambient  air 

Qc,i 

heat  rejected  by  the  i-th  subsection  of  TE  mod¬ 

c 

cold  junction  of  TE  module 

ule  (W) 

coolant 

coolant 

Qconv 

convective  heat  loss  (W) 

F 

Fresnel  lens 

Qcool 

heat  removed  by  the  cooling  system  (W) 

h 

hot  junction  of  TE  module 

Qh 

heat  absorbed  by  the  hot  junction  of  the  TE 

HC 

heat  collector 

module  (W) 

HS 

heat  sink 

Qhj 

heat  absorbed  by  the  i-th  subsection  of  TE  mod¬ 

i 

the  i-th  subsection  of  TE  module 

ule  (W) 

PV 

photovoltaic  solar  cell 

Qrad 

radiative  heat  loss  (W) 

PV-TE 

PV-TE  hybrid  system 

R 

electrical  resistance  (Q) 

TE 

TEG  subsystem 

Rl 

external  load  (Q) 

ks 

split  at  cutoff  wavelength,  ks 

Rs 

series  resistance  of  solar  cell  (Q) 

Rte 

resistance  of  TE  module  (Q) 

The  spectrum  response  of  many  solar  cells  is  most  effi¬ 
cient  when  photon  energy  is  close  to  solar  cell  band-gap 
energy.  Photons  of  energy  smaller  or  larger  than  band- 
gap  energy  can  only  be  partly  utilized.  In  order  to  make  full 
utilization  of  solar  energy  in  the  broad  solar  spectrum 
wavelength  range,  the  notion  of  spectrum  splitting  PV- 
TE  hybrid  systems  is  developed.  By  combining  the  photo¬ 
voltaic  module  with  the  solar  thermoelectric  generator 


(TEG),  photons  out  of  range  of  a  given  solar  cell’s  narrow 
absorption  wavelength  range  can  be  directed  to  TEG 
modules,  generating  electricity  through  thermoelectric 
effect.  The  energy  conversion  efficiency  would  benefit  from 
this  procedure,  and  the  heat  dissipated  by  the  PV  module  is 
reduced.  So  far,  several  researches  have  been  done  on  spec¬ 
trum  splitting  PV-TE  hybrid  systems.  Zhang  et  al.  (2005) 
designed  a  PV-TE  power  generating  system  for  a  project 
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of  “Nano  and  graded  thermoelectric  materials/photovol¬ 
taic-thermoelectric-wind  power  generation” .  V orobiev 
et  al.  (2006a,  2006b)  discussed  the  possibility  of  spectrum 
splitting  PV-TE  hybrid  systems,  showing  that  this  hybrid 
system  could  be  practical  and  efficient.  Kraemer  et  al. 
(2008)  developed  a  general  optimization  methodology  for 
spectrum  splitting  PV-TE  hybrid  systems.  In  his  study, 
an  optimum  methodology  for  cutoff  wavelength  of  such 
hybrid  systems  was  developed  and  the  power  fractions  of 
each  segment  were  analyzed.  It  was  revealed  that  the  short 
wavelength  region  which  should  be  directed  to  the  solar 
TEG  is  usually  only  a  small  portion  of  the  total  solar  spec¬ 
trum  and  could  be  neglected.  Fleurial  (2009)  also  men¬ 
tioned  that  the  PV-TE  hybrid  solar  power  system  based 
on  splitting  the  solar  energy  spectrum  could  help  maximize 
the  conversion  efficiency  and  improve  thermal  management 
of  highly  concentrated  solar  systems. 

However,  the  influence  of  concentration  ratio  and  cooling 
condition,  which  is  significant  for  such  systems  working 
under  high  concentrations,  was  not  discussed  in  previous 
studies.  For  practical  application  of  the  PV-TE  hybrid  sys¬ 
tem,  we  developed  a  numerical  model  in  this  paper,  taking 
into  consideration  the  solar  energy  concentration  and  the 
semiconductor  material  working  temperature.  This  model 
is  helpful  to  study  the  influences  and  relations  of  these 
parameters.  For  instance,  a  GaAs-CoSb3  hybrid  PV-TE  sys¬ 
tem  was  analyzed  using  this  model,  relationships  between 
cutoff-wavelength,  concentration  ratio  and  heat  sink  perfor¬ 
mance  were  studied,  and  guidelines  for  the  PV-TE  hybrid 
system  design  and  optimization  were  provided. 

2.  Spectrum  splitting  concentration  PV-TE  hybrid  system 

Spectrum  splitting  concentration  photovoltaic-thermo¬ 
electric  (PV-TE)  hybrid  systems  consist  of  concentrators, 
photovoltaic  solar  cells,  thermoelectric  generators,  spectral 
beam  splitters  and  cooling  systems.  As  described  by  Imenes 
and  Mills  (2004),  many  different  types  of  spectrum  splitting 
techniques  have  been  proposed  over  the  years  for  solar  sys¬ 
tems.  (1)  Dichroic  splitters:  dichroic  spectral  beam  splitters 
are  mirrors  or  prism  assemblies  with  dielectric/metal-dielec- 
tric  optical  coatings  (Borden  et  al.,  1982;  Barnett  et  al., 
2007).  Depending  on  the  coating  being  used,  reflection/ 
transmission  ratios  may  differ  in  function  of  the  wavelength. 
(2)  Prisms:  a  prism  will  break  incident  light  up  to  its  constit¬ 
uent  spectral  beams.  It  has  been  applied  in  spectral  beam 
split  systems  for  different  band-gap  solar  cells  (Converse, 
1997).  (3)  Absorption  filters:  it  typically  consists  of  liquid 
contained  glass  plates  or  tubes.  The  liquid  to  which  various 
inorganic  or  organic  compounds  have  been  added  absorbs 
wavelengths  of  light  above  the  cell’s  band-gap  and  converts 
the  energy  to  heat,  which  could  be  utilized  by  other  methods 
(Hamdy  and  El-Hefnawi,  1990;  Sabry  et  al.,  2002).  Other 
possible  solutions  to  split  sun  light  by  spectrum  wavelengths 
including  luminescent  concentrators  (Galluzzi  and  Scafe, 
1984;  Van  Sark  et  al.,  2008;  Fisher  and  Biddle,  2011)  and 
holographic  concentrators  (Ludman  et  al.,  1997;  Hull 


et  al.,  1987),  these  devices  that  both  concentrate  and  spec¬ 
trally  split  solar  light  can  be  used  in  low-concentration 
non-tracking  solar  energy  systems. 

For  the  spectrum- splitting  hybrid  system  with  separated 
PV  receiver  and  TE  generator,  it’s  possible  to  employ  all 
these  types  of  splitters  by  using  different  system  configura¬ 
tions.  Taking  a  PV-TE  hybrid  system  with  a  Fresnel  lens 
and  a  dichroic  mirror  for  example  (Fig.  la),  incident  sunlight 
is  firstly  concentrated  by  the  Fresnel  lens.  The  concentrated 
sunlight  is  focused  on  the  dichroic  beam  splitter  and  split 
into  two  segments  at  cutoff-wavelength  ks.  Irradiance  of 
spectrum  wavelength  below  As,  in  which  the  concentrated 
photovoltaic  (CPV)  cell  is  of  good  performance,  is  directed 
to  CPV  cells,  while  the  other  part  of  irradiance  is  directed 
to  TE  generators.  Because  solar  cells  working  under  high 
radiation  flux  suffer  from  high  temperatures,  which  are 
harmful  to  cell  efficiency  and  operational  life  span,  a  cooling 
system  is  used  to  remove  the  dissipated  heat  from  the  CPV 
subsystem.  For  the  TEG  subsystem,  the  dissipated  heat  is 
also  removed  by  cooling  systems  which  can  provide  a  con¬ 
stant  temperature  difference  for  the  thermoelectric  effect  of 
the  solar  TEG.  The  segmentation  of  solar  spectrum  is  shown 
in  Fig.  2.  As  cutoff  wavelength  As  increases,  power  fractions 
of  CPV  subsystem  and  TEG  subsystem  show  opposite 
trends.  In  this  paper,  we  use  the  AM1.5D  solar  spectrum, 
which  is  defined  by  ASTM  G  173-03  (2003)  for  solar  concen¬ 
trator  work. 

3.  Numerical  module 

The  spectrum  splitting  PV-TE  system  normally  works 
under  high  solar  concentration  and  high  temperature  con¬ 
ditions.  In  order  to  determine  the  thermal  and  electrical 
performance  of  the  hybrid  system,  a  1 -dimensional  (ID) 
numerical  model  was  established  for  simulation.  In  this 
model,  we  assumed  that:  (1)  the  solar  cell  and  heat  collector 
are  of  the  same  effective  area,  APV  =  AHC ;  (2)  the  thermal 
insulation  materials  of  the  thermoelectric  generator  are  of 
very  high  performance  and  the  heat  transfer  in  these  mate¬ 
rials  is  neglected;  (3)  the  same  densely  packaged  active 
cooling  systems  are  adopted  by  both  the  CPV  and  the 
TEG;  (4)  As  solar  cells  are  very  thin  and  concentrating 
solar  cells  usually  have  high  thermal  conductivity,  temper¬ 
ature  gradients  of  solar  cells  are  neglected  in  numerical 
modeling  for  simplification. 

Fig.  lb  shows  the  heat  transfer  and  energy  conversion  of 
the  PV-TE  hybrid  system.  The  solar  irradiance  Ein  is  con¬ 
centrated  and  split  into  two  parts  at  cutoff  wavelength  ks. 
The  energy  directed  to  the  TEG  subsystem  and  the  CPV 
subsystem  can  be  calculated  with  the  following  equations, 


4000nm 


ETe  —  /  rjoptAHCCgF(X)d/ 1 

Jls 

(i) 

rls 

Epv  =  /  r\optAPVCgF(X)ds l 

(2) 
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Fig.  1.  The  PV-TE  hybrid  system:  (a)  Schematic  diagram  of  the  hybrid  system;  (b)  Energy  flow  chart  of  the  hybrid  system. 


where  Ete  is  the  solar  energy  directed  to  the  TEG  subsys¬ 
tem  (from  ks  to  4000  nm);  EPV  is  the  solar  energy  directed 
to  the  CPV  subsystem  (from  280  nm  to  A5);  F( A)  is  the  spec¬ 
tral  irradiance  in  AM1.5D  condition,  and  rjopt  is  the  effi¬ 
ciency  of  the  optical  system.  Geometrical  concentration 
ratio  Cg  is  determined  by  ratio  of  the  Fresnel  lens  area 
and  the  device  surface  area.  As  mentioned  above,  it  is 
equivalent  for  both  TEG  and  CPV:  Cg  =  AHCI AF  —  APVI 
A F.  Ahc,  Apv  and  AF  are  the  surface  area  of  the  heat  col¬ 
lector,  the  solar  cell  and  the  Fresnel  lens  respectively. 

For  the  TEG  subsystem,  most  of  the  radiant  energy  Ete 
is  absorbed  by  the  heat  collector  and  converted  into  ther¬ 
mal  energy.  The  heat  losses  of  the  heat  collector  are  caused 
by  radiation  and  natural  convection.  The  radiative  heat 
loss  of  the  heat  collector  is, 


Qrad,TE  ~  ZhcAHcCsb(TAhc  ~  T\ir)  (3) 

where  THC  is  the  surface  temperature  of  the  heat  collector, 
T air  is  the  ambient  air  temperature,  sHC  is  the  emissivity  of 
the  heat  collector  surface  coating,  gSb  is  the  Stefan-Boltz- 
mann  constant. 

And  the  heat  loss  through  convection  of  the  heat  collec¬ 
tor  surface  is, 

Qconv.TE  =  KAhc{THC  ~  T air)  (4) 

where  hnc  is  the  natural  convection  heat  transfer  coefficient. 

In  the  case  of  single-phase  forced  convection  cooling, 
the  heat  removed  by  the  cooling  system,  which  is  also  the 
heat  rejected  by  the  cold  junction  of  the  thermoelectric 
(TE)  module,  can  be  obtained  by, 

Q  cool, IE  =  hCOolAHs(T c  —  T coolant  )  (5) 

where  Tc  and  Tcooiant  are  the  mean  cold-junction  and 
coolant  temperatures,  hcooi  is  the  heat  transfer  coefficient 


Fig.  2.  Segmentation  of  AM1.5D  spectrum  in  two  regions  for  the  CPV 
subsystem  and  TEG  subsystem.  The  inset  shows  the  converted  power 
fraction  of  each  subsystem  when  cutoff  wavelength  point  increases. 

which  is  commonly  used  for  comparing  the  heat  transfer 
characteristics  of  cooling  systems,  A HS  is  the  effective  sur¬ 
face  area  for  heat  transfer  between  the  heat  sink  and  the 
cold-junction  of  TE  module,  AHS-  A HC. 

The  heat  absorbed  by  the  hot  junction  of  the  TE  module 
can  be  expressed  as, 

Qh,TE  =  UhcEtE  ~  Qrad,HC  ~  Qconv,HC  (6) 

where  ocHC  is  the  absorptivity  of  the  heat  collector  surface 
coating. 

Hence,  the  generated  power  of  the  TE  module  is, 

Ete  —  Qh,TE  ~  Q cool, te  (7) 
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For  the  CPV  subsystem,  the  heat  transferred  by  radia¬ 
tion  and  convection  on  solar  cell  surface,  and  the  heat 
removed  by  cooling  system  are  expressed  in  similar  formu¬ 
las  as  the  TEG  subsystem, 

Qrad,PV  =  ZpvApVGSB(TpV  —  TAIR)  (8) 

Qconv,PV  =  hncApv(TpV  ~  T AIR)  (9) 

Qcool,PV  =  hCoolA-Hs{T PV  —  T coolant  )  (10) 

where  TPVis  the  temperature  of  the  solar  cell  and  &PV\s  the 
emissivity  of  the  solar  cell  surface. 

Based  on  the  energy  balance  analysis,  the  generated 
power  of  the  CPV  subsystem  can  be  expressed  as, 

PpV  =  dpvEpv  —  Qradfv  ~  QCOnvPV  ~  QcoolPV  (H) 

where  aPV  is  the  absorptivity  of  the  solar  cell  surface. 

The  total  generated  power  of  the  PV-TE  hybrid  system 
is  the  combination  of  the  CPV  subsystem  output  power 
and  the  TEG  subsystem  output  power, 

P pv-TE  =  P pv  A-  Pte  (12) 

And  the  conversion  efficiency  of  the  PV-TE  hybrid  sys¬ 
tem  is, 

Vpv-te  =  P pv-te / {Af  •  Ein)  (13) 

Eqs.  (1-13)  describe  the  energy  conversion  of  the  PV- 
TE  hybrid  system.  Base  on  the  fundamental  equations 
described  above,  models  for  other  layouts  can  be  easily 
extended.  However,  the  temperature  distribution  of  the 
hybrid  system  is  still  unknown,  it’s  necessary  to  establish 
the  thermal-electric  models  for  the  solar  cell  and  the  TE 
module  separately  and  solve  the  models  with  Eqs.  (1-13) 
simultaneously.  As  the  performance  of  thermoelectric 
material  will  be  highly  affected  by  the  temperature  gradient 
under  high  concentrations,  we  use  a  finite  element  model  to 
simulate  the  TE  generator  in  this  paper  (Li  et  al.,  2010). 
And  for  the  solar  cell,  a  temperature  corrected  Cheknane 
model  (Cheknane  et  al.,  2006)  is  introduced  for  a  better 
prediction  of  solar  cell  efficiency.  Details  and  explanations 
of  how  we  set  up  these  models  of  subsystems  are  described 
in  Appendix. 

Based  on  the  algorithm  above,  a  calculation  program 
has  been  developed  for  the  PV-TE  hybrid  system.  As 
shown  in  Fig.  3,  solar  radiation  is  split  at  the  cutoff  wave¬ 
length  and  introduced  to  the  CPV  and  the  TEG  subsystems 
as  energy  input.  For  the  TEG  subsystem,  an  iterative  algo¬ 
rithm  has  been  developed.  Providing  initial  values  of  tem¬ 
perature  distribution,  we  can  calculate  the  temperature 
dependent  material  properties  of  each  finite  element:  elec¬ 
trical  resistance  Rh  thermal  conductivity  Kt  and  Seeback 
coefficient  a*.  These  parameters  are  substituted  into  energy 
balance  equations  of  the  TEG  subsystem  to  obtain  a  new 
temperature  distribution  by  using  discrete  Newton  algo¬ 
rithm.  The  calculated  new  temperatures  are  used  for  the 
next  iteration  until  the  values  of  temperature  distribution 
converge.  For  the  CPV  subsystem,  we  deduced  the  solar 
cell  electrical  performance  under  1  sun  and  split- spectrum 


Fig.  3.  Flow  chart  of  calculation  procedure  for  the  hybrid  system. 


conditions  from  the  solar  cell  external  quantum  efficiency 
(EQE)  and  solar  cell  open-circuit  voltage  Voc.  These 
parameters  of  solar  cell  performance  are  substituted  to 
the  temperature  corrected  Cheknane  model  to  predict  elec¬ 
trical  performance  of  the  solar  cell  under  concentration 
ratio  Cg.  Based  on  the  energy  balance  equation  of  the  solar 
cell,  working  temperature  TPV  can  be  obtained.  After  that, 
the  temperature  distribution  of  the  hybrid  system  is 
revealed,  output  power,  conversion  efficiency  and  other 
electrical  and  thermal  results  of  the  TEG  and  CPV  systems 
can  be  obtained.  Finally,  we  combine  the  results  calculated 
from  the  TEG  and  the  CPV  subsystems,  and  the  electrical 
and  thermal  performance  of  the  hybrid  system  can  be  fig¬ 
ured  out. 


4.  Results  and  discussions 

In  order  to  determine  the  influence  of  cutoff  wavelength, 
concentration  ratio  and  heat  dissipation  on  the  PV-TE 
hybrid  system,  a  PV-TE  hybrid  system  with  a  single  junc¬ 
tion  GaAs  solar  cell  (band  gap  Eg  =  1 .43  eV)  and  a  high 
performance  skutterudites  (figure  of  merit  ZT  —  1.4  at 
800  K)  thermoelectric  module  is  theoretically  simulated. 
Performance  and  material  properties  of  the  solar  cell  and 
thermoelectric  semiconductor  are  obtained  from  literatures 
of  Algora  et  al.  (2001)  and  Li  et  al.  (2009),  both  n-type  and 
p-type  legs  of  TE  module  are  assumed  to  be  the  same  ther¬ 
moelectric  properties.  Other  parameters  using  in  this  model 
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Table  1 


Parameters  used  in  the  PV-TE  hybrid  model. 


Symbol 

Description 

Value 

Symbol 

Description 

Value 

Optical  system 

Vopt 

optical  efficiency 

85% 

GaAs  Solar  cell 

nf 

Ideality  factor 

1.9 

ft, 

Temperature  coefficient 

-0.002  K  1 

ttpy 

Absorptivity 

0.95 

fisc 

Temperature  coefficient 

-0.0007  K-1 

Spy 

Apy 

Emissivity 

Surface  area 

0.85 

1  x  1  mm2 

Pvoc 

Temperature  coefficient 

-0.002  K“‘ 

Heat  collector 

Khc 

Thermal  conductivity 

0.2  W/K 

rjHC 

Absorptivity 

0.9 

Arc 

Surface  area 

1  x  1mm2 

£hc 

Emissivity 

0.08 

CoSb3  TE 
module 

An 

Ap 

n-type  leg  cross-section 

area 

p-type  leg  cross-section 

area 

0.4  x  0.4mm2 

0.4  x  0.4  mm2 

L 

n-type  and  p-type  leg  length 

3.6  mm 

Others 

Tair 

Ambient  air  temperature 

25  °C 

'A coolant 

Coolant  temperature 

25  °C 

hP 

Planck  constant 

6.626068  x  10“34m2kg/s 

C 

Speed  of  light  in  a  vacuum 

2.9979  x  108  m/s 
1.602176  x  10“19 

O SB 

Stefan-Boltzmann 

constant 

5.67  x  10“8  W/(m2  K4) 

Q 

Elementary  charge 

1.602176  x  10~19C 

kp 

Boltzmann  constant 

1.38065  x  10"23m2kg/ 

(s2  K) 

hnc 

Natural  heat  transfer 
coefficient 

5  W/(m2  K) 

are  given  in  Table.  1.  All  calculations  are  limited  by  the 
assumed  operating  temperature,  of  which  the  maximum 
operating  temperature  of  CoSb3  TE  module  is  800  K  and 
that  of  the  GaAs  solar  cell  is  500  K. 

4.1.  Influence  of  cutoff  wavelength 

The  generated  power  of  the  PV-TE  hybrid  system  when 
cutoff-wavelengths  are  ranged  from  300-1800  nm  is  shown 
in  Fig.  4,  under  the  concentration  ratio  of  200  suns  and  heat 
transfer  coefficient  of  10,000  W/m2  K  condition.  Since  the 
energy  partition  is  determined  by  cutoff-wavelength,  output 
power  of  the  hybrid  system  and  subsystems,  P pv-te,  Ppv 
and  PTE,  is  obviously  varied  with  ks.  The  maximum  hybrid 
system  output  power  is  achieved  at  cutoff  wavelength 
850-950  nm  in  this  condition.  Output  power  of  the  TEG 


Cutoff  Wavelength  A  (nm) 

Fig.  5.  Working  temperature  of  the  hybrid  system  vs.  cutoff  wavelength, 
Cg  =  200,  h  =  10,000  W/m2  K_1. 


Fig.  4.  Generated  power  of  the  hybrid  system  vs.  cutoff  wavelength, 
Cg  =  200,  h=\  0,000  W/m2  K“ 1 . 


subsystem  reduces  when  cutoff  wavelength  rises,  it  falls  dra¬ 
matically  in  ks  =  400-800  nm  because  of  the  energy  distri¬ 
bution  in  the  solar  spectrum.  For  the  CPV  subsystem, 
output  power  PPV increases  rapidly  until  it  reaches  the  max¬ 
imum  point,  then  reduces  slightly.  Investigation  of  the  sub¬ 
system’s  working  temperature  (Fig.  5)  shows  that  the  solar 
cell  temperature  TPV  increases  slowly  with  ks.  And  the  tem¬ 
perature  difference  between  the  hot-  and  cold-junction  of 
the  TE  module,  A  T=  Th  —  Tc,  decreases  rapidly  with  2, 
which  leads  to  the  TEG  output  power  reduction. 

For  a  further  study  of  the  relationship  between  concen¬ 
tration  ratios  and  optimized  cutoff-wavelengths,  Cg  200- 
600  were  calculated  using  the  same  process.  As  can  be  seen 
in  Fig.  6,  detailed  results  have  shown  that  when  the  system 
efficiency  increases  with  concentration  ratio,  the  optimized 
cutoff  wavelength  keeps  steady.  This  is  because  EQE  of  the 
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Fig.  6.  Hybrid  system  efficiency  vs.  cutoff  wavelength  for  different 
concentration  ratio,  h  =  10,000  W/m2  K-1. 


GaAs  solar  cell  falls  sharply  when  energy  of  photons  is  less 
than  band  gap  Eg  (photos  of  wavelength  larger  than  hP-d 
Eg  =  867  nm  interact  only  weakly  with  the  solar  cell). 
Consequently,  we  accept  X  =  900  nm  as  the  optimized  cut- 
off-wavelength  of  this  spectrum  splitting  PV-TE  hybrid 
system.  At  this  point,  incident  direct  solar  energy 
Ein  =  900  W/m2  is  divided  into  596.52  W/m2  for  the  CPV 
subsystem  and  303.62  W/m2  for  the  TEG  subsystem,  which 
means  that  the  incident  energy  partition  for  the  CPV  sub¬ 
system  is  66.3%,  and  for  the  TEG  subsystem  is  33.7%. 

Generated  power  of  the  hybrid  system  working  under 
cutoff-wavelength  Xs  =  900  nm  is  highly  dependent  on  the 
CPV  subsystem  because  this  subsystem  converts  the  bigger 
portion  of  solar  energy  and  has  a  higher  conversion  effi¬ 
ciency.  As  calculated  in  Fig.  7,  the  hybrid  system  generates 
power  Pout  and  the  TEG  subsystem  generates  power  PTE 
increasing  with  the  concentration  ratio.  The  partition  of 
TEG  subsystem  generated  power,  Pte/Ppv-te ,  also 
increases  with  Cg ,  but  it’s  still  less  than  12%  at  800  suns. 
Generated  power  of  the  hybrid  system  is  mostly  provided 
by  the  GaAs  solar  cell. 


Fig.  7.  Generated  power  and  TEG  subsystem  contribution  vs.  concen¬ 
tration  ratio,  h  =  10,000  W/m2  K-1. 


4.2.  Influence  of  concentration  ratio 

When  the  PV-TE  hybrid  system  works  at  the  optimized 
cutoff  wavelength  Xs  =  900  nm,  electrical  and  thermal  per¬ 
formance  of  the  hybrid  system  will  significantly  change 
with  the  concentration  ratio  of  the  optical  system  and  the 
performance  of  the  heat  sink.  In  this  study,  we  adopted 
water  (the  typical  heat  transfer  coefficient  for  water  is  ran¬ 
ged  from  500  to  10,000  W/m2  K)  as  the  coolant  of  the  heat 
sink  and  several  different  values  of  heat  transfer  coefficient 
were  considered  in  the  calculation. 

Fig.  8  shows  calculated  temperatures  of  the  hybrid  sys¬ 
tem  according  to  the  concentration  ratio  for  several  values 
of  hcooi.  The  solar  cell  temperature  TPV ,  the  hot-  and  cold- 
junction  temperatures  of  TE  module,  Th  and  Tc,  and  their 
temperature  difference  A  T  increase  linearly  with  concentra¬ 
tion  ratio  Cg.  It  can  be  observed  from  Fig.  8a  and  b  that  the 
maximum  concentration  ratios  are  limited  by  the  solar  cell 
working  temperature  for  hcooi  =  500-1000  W/m2  K.  For 
hcooi  =  2000  to  10,000  W/m2  K  the  hot-junction  tempera¬ 
ture  of  the  TE  module  is  the  limiting  factor.  Curves  of 
Fig.  8d  are  almost  at  the  same  slope  for  all  hcooi  values, 
which  means  that  the  hot-  and  cold-junction  temperature 
difference  A  T  is  almost  independent  of  hcooi •  Additionally, 
TPV ,  Th  and  Tc  of  a  hybrid  system  with  high  heat  transfer 
coefficient  hcooi  are  increasing  slower  than  that  of  a  small 
hCooi •  F or  the  hybrid  system  with  a  larger  heat  transfer  coef¬ 
ficient,  higher  concentration  ratios  could  be  achieved. 

Fig.  9  shows  the  efficiency-C^  curves  for  different  cool¬ 
ing  conditions.  As  shown  in  Fig.  9a,  for  cooling  conditions 
of  hcooi  =  500,  1000  and  2000  W/m2  K,  the  hybrid  system 
efficiency  rj  pv-te  increases  according  to  concentration  ratio 
Cg ,  reaching  its  maximum  point  within  concentration  limits 
and  then  falling  down.  For  heat  transfer  coefficient 
hcooi  =  5000,  10,000  and  oc  W/m2  K,  tjfv-te  keeps  increas¬ 
ing  with  concentration  ratio.  Fig.  9b  shows  the  CPV  sub¬ 
system  efficiency  according  to  the  concentration  ratio. 
Similar  to  the  first  type  of  the  hybrid  system  efficiency 
curves,  the  CPV  subsystem  efficiency  tjpv  exhibits  a  maxi¬ 
mum  with  respect  to  the  concentration  ratio.  The  reduction 
of  the  CPV  subsystem  efficiency  rjpv  at  high  concentration 
is  due  to  the  rising  temperature  and  the  series  resistance  of 
the  solar  cell.  The  solar  cell  efficiency  in  the  PV-TE  hybrid 
system  is  much  higher  than  in  the  solar  cell  only  system. 
Under  hcooi  =  °o  condition,  maximum  solar  cell  efficiency 
of  the  PV-TE  hybrid  system  is  about  40%,  while  it  is  only 
26%  for  GaAs  solar  cell.  As  shown  in  Fig.  9c,  the  TE  sub- 

o 

system  efficiency  rj TE  increases  with  Cg.  Since  the  TEG  sub¬ 
system  efficiency  is  mainly  determined  by  the  temperature 
difference  between  hot-  and  cold-junctions,  curves  of  r]TE 
are  independent  of  hcooi ,  they  also  show  a  linearly  increas¬ 
ing  trend  with  Cg. 

4.3.  Influence  of  heat  transfer  coefficient 

As  analyzed  before,  the  PV-TE  hybrid  system  efficiency 
will  increase  with  the  heat  transfer  coefficient  of  the  heat 
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Fig.  8.  Temperature  vs.  concentration  ratio  for  different  heat  transfer  coefficient,  X  =  900  nm.  (a)  Solar  cell  temperature;  (b)  TE  hot-junction  temperature; 
(c)  TE  cold-junction  temperature;  (d)  AT  of  TE  hot-  and  cold-  junctions. 


sink,  but  higher  performance  heat  sinks  usually  need  bigger 
pumps  and  more  complex  designs,  which  will  lead  to  high 
parasite  energy  and  expensive  costs.  Therefore,  the  influ¬ 
ence  of  the  heat  transfer  coefficient  on  the  hybrid  system 
efficiency  is  also  analyzed  in  this  study.  As  shown  in 
Fig.  10,  for  concentration  ratio  Cg  =  200,  400,  600,  800 
and  cutoff  wavelength  Xs  =  900  nm,  when  the  heat  transfer 
coefficient  hcooi  increases  from  500  to  3000  W/m2  K,  system 
efficiency  of  the  hybrid  system  fj  pv-te  increases  rapidly,  but 
for  the  heat  transfer  coefficient  hcooi  larger  than  3000  W/ 
m2  K,  Opv-te  will  not  change  much  with  increasing  hcooi. 

Fig.  1 1  shows  how  the  hybrid  system  efficiency  and  the 
concentration  ratio  vary  with  the  heat  transfer  coefficient 
when  the  system  achieves  its  maximum  efficiency.  A  PV- 
only  system  which  uses  the  same  GaAs  solar  cell  is  also  cal¬ 
culated.  The  maximum  conversion  efficiency  of  the  PV-TE 
hybrid  system  and  the  PV-only  system  increase  with  the 
heat  transfer  coefficient,  but  the  rate  of  increase  slows 
down.  Here  we  use  Cg\n  to  represent  the  concentration 


ratio  when  the  system  achieves  its  maximum  efficiency  at 
a  certain  heat  transfer  coefficient.  This  figure  also  shows 
that  Cg \nMAX  increases  almost  linearly  with  heat  transfer 
coefficient  in  a  certain  range  (For  the  PV-TE  hybrid 
system,  hcooi  ^4500  W/m2  K;  For  the  PV-only  system, 
hcooi  5H0000  W/m2  K).  In  this  paper,  we  define  the  slope 

of  Cg\riMAx-hcooi  curves,  y  =  dC^x,  as  the  feasibility  of  sys¬ 
tem  working  under  concentration.  The  PV-TE  hybrid  sys¬ 
tem  of  y pv-te  ~  0.193  is  more  suitable  for  working  under 
high  concentrations  than  the  PV-only  system  of  ypv-oniy  ~ 
0.027.  In  the  range  of  hcooi  >  4500  W/m2-K.  the  concentra¬ 
tion  ratio  of  the  PV-TE  hybrid  system  is  restricted  by  the 
hot-junction  temperature  of  the  TE  module,  the  maximum 
efficiency  point  being  achieved  at  Th  =  800  K. 

The  thermal  management  of  the  hybrid  system  is  a  sig¬ 
nificant  challenge.  For  densely  packaged  PV-TE  hybrid 
systems  working  under  high  concentration  ratio  (>300 
suns),  it’s  not  a  good  solution  to  use  regular  cooling 


X.  Ju  et  al.  I  Solar  Energy  86  (2012)  1941-1954 


1949 


(a)  (b) 


(c) 


h-  co 

-9- 

h-500W/m2K 

— B- 

h=1000W/m2K 

— e- 

h=2000W/m2K 

-  h=5000W/m2K 

— A— 

h=10000W/m2K 

Fig.  9.  Efficiency  vs.  concentration  ratio,  X  =  900  nm.  (a)  Hybrid  system  efficiency;  (b)  CPV  subsystem  efficiency;  (c)  TEG  subsystem  efficiency. 


Fig.  10.  Hybrid  system  efficiency  vs.  heat  transfer  coefficient  of  heat  sink, 
X  =  900  nm. 


Fig.  11.  Concentration  ratio  and  system  efficiency  vs.  heat  transfer 
coefficient,  when  the  PV-TE/PV-only  system  efficiency  reaches  the 
maximum  point. 
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Table  2 


Performance  of  the  PV-TE  hybrid  system  in  optimized  range. 


h  (W/m2  KT1) 

Cg 

Pout  (W) 

Ppv(  W) 

Pte  (W) 

Pte/Pou,  (%) 

Tpv(  K) 

Th  (K) 

TA  K) 

AT  (K) 

V] PV-TE  (%) 

3000 

550 

0.1318 

0.1202 

0.0115 

8.73 

352.2 

694.0 

336.1 

357.9 

26.62 

4500 

770 

0.1905 

0.1697 

0.0208 

10.92 

348.4 

804.7 

332.3 

468.5 

27.49 

Fig.  12.  Comparison  of  the  efficiency  between  the  PV-only  system  and  the 
PV-TE  hybrid  system. 


method.  Although  the  calculated  heat  transfer  coefficient  is 
between  3000  and  4500  W/m2  K,  we  should  take  into  con¬ 
sideration  the  system  working  at  open  circuit  and  DNI 
(Direct  Normal  Irradiance)  =  1300  W/m2  condition.  For 
the  safety  and  ruggedness  of  the  PV-TE  hybrid  system 
working  under  extreme  terrestrial  conditions,  the  cooling 
system  of  hcooi  >  10000  W/m2  K  is  necessary.  According 
to  the  reported  values  reviewed  by  Royne  et  al.  (2005), 
potential  cooling  methods  include  microchannel  cooling, 
jet  impingement  cooling  and  phase  change  cooling.  How¬ 
ever,  it  should  be  noticed  that  the  phase  change  cooling 
methods  usually  use  toxic  coolants,  which  will  lead  to  addi¬ 
tional  complex  design  and  expensive  costs. 

4.4.  Optimization  of  the  PV-TE  hybrid  system 

Guidelines  for  design  and  optimization  are  summarized 
from  the  correlations  developed  in  this  paper.  For  the  spec¬ 
trum  splitting  PV-TE  hybrid  system  design,  the  optimal 
device  design  procedure  can  be  carried  out  in  the  following 
steps: 

(1)  Determine  the  type  of  solar  cells  and  TE  generators. 
Selected  materials  should  be  suitable  for  working 
under  concentrations.  Under  high  concentrated  sun¬ 
light  (>300  suns),  III-V  based  solar  cells  can  be  used 
to  produce  electricity  in  the  CPV  subsystem,  while 
skutterudites,  TAGS,  PbTe  or  LAST  alloy  thermo¬ 
electric  modules  are  potential  choices  of  the  TEG 
subsystem. 


T  coolant 


Rl 

Load 


Fig.  13.  Schematic  representation  of  a  thermal  module  of  a  single-couple 
TEG. 


(2)  Find  the  optimized  cutoff-wavelength  ks.  It’s  usually 
decided  by  the  band  gap  Eg  of  solar  cell,  ks^hP-cl 
Eg.  Spectral  beam  splitter  could  be  selected  according 
to  the  cutoff  wavelength.  Dichroic  mirrors  and  prisms 
are  good  choices,  since  they  will  make  the  system 
design  more  simple. 

(3)  Determine  the  cooling  system  of  the  hybrid  system.  It  is 
better  to  use  a  cooling  system  of  which  the  heat  transfer 
coefficient  is  1-2  times  higher  than  the  optimized  range 
of  hcooi  in  practical  application.  Usually,  cooling  sys¬ 
tems  of  the  heat  transfer  coefficient  >  10000  W/m2  K 
is  necessary  for  such  hybrid  system  working  under  high 
concentrations.  The  cooling  system  should  be  light  in 
weight,  since  the  system  usually  tracking  the  sun. 

(4)  Figure  out  the  concentration  ratio  when  the  system 
achieves  its  maximum  efficiency  at  the  optimized 
range  of  hcooi.  The  concentration  system  design  of 
the  hybrid  system  is  determined  by  the  concentration 
ratio,  both  point-focus  or  line-focus  solar  collectors 
can  be  used  to  concentrate  sunlight. 

(5)  Determine  the  other  system  configurations  including 
the  tracking  system  and  power  control  system. 

For  the  GaAs-CoSb3  hybrid  system  analyzed  in  this 
work,  the  cutoff  wavelength  ks  =  900  nm.  As  shown  in 
Table. 2,  the  optimized  heat  transfer  coefficients  are  ranged 
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from  3000  to  4500  W/m1 2 3  K  and  the  optimized  concentra¬ 
tion  ratio  of  the  hybrid  system  is  from  550  to  770  suns. 
In  this  optimized  working  region,  power  fraction  of  TEG 
subsystem  is  ranged  from  8.73%  to  10.92%,  and  the  hybrid 
system  efficiency  is  from  26.62%  to  27.49%.  The  solar  cell 
temperature  keeps  at  around  350  K,  which  is  the  proper 
working  temperature  for  GaAs  solar  cells.  The  temperature 
difference  between  the  hot-  and  cold-junction  of  the  TE 
module  is  about  350-470  K.  Fig.  12  shows  the  comparison 
of  system  efficiency  between  the  PV-TE  hybrid  system  and 
the  PV-only  system.  When  systems  work  in  the  maximum 
efficiency  region,  PV-TE  hybrid  systems  could  achieve 
higher  efficiency  than  PV-only  systems.  What  is  more 
important,  the  concentration  ratio  of  the  PV-only  system 
is  only  110-160  suns  in  the  optimized  region.  It  is  much 
lower  than  PV-TE  hybrid  system. 

It  can  also  be  deduced  from  Table. 2  and  Fig.  9a  that 
there  are  two  basic  ways  to  improve  the  hybrid  system  effi¬ 
ciency.  One  is  to  increase  the  concentration  ratio  of  the  sys¬ 
tem.  Assuming  that  the  heat  dissipation  rate  of  the  heat 
sink  is  large  enough,  the  higher  the  concentration  ratio, 
the  greater  the  hybrid  system  efficiency  can  be  achieved. 
The  contribution  of  TEG  subsystem  will  increase  with 
the  concentration  ratio.  The  other  is  to  improve  the  perfor¬ 
mance  of  the  cooling  system.  This  method  helps  to  reduce 
the  temperature  of  the  solar  cell,  in  other  words,  helps  to 
reduce  the  efficiency  reduction  caused  by  increasing  tem¬ 
perature.  The  potential  direction  of  system  development 
is  to  use  high  temperature  thermoelectric  modules,  elevat¬ 
ing  the  concentration  ratio  and  simultaneously  choosing 
a  proper  cooling  system. 

5.  Conclusion 

Energy  conversion  and  heat  transfer  process  of  the  spec¬ 
trum  splitting  concentration  photovoltaic-thermoelectric 
hybrid  system  is  investigated  in  this  paper,  and  an  energy- 
based  numerical  model  for  PV-TE  hybrid  systems  is 
presented.  This  model  is  used  to  analyze  the  thermal  and 
electrical  performance  of  hybrid  systems  with  GaAs  solar 
cells  and  skutterudites  CoSb3  solar  TE  generators.  Results 
show  that  the  spectrum  splitting  PV-TE  hybrid  system 
has  several  potential  advantages  over  a  PV-only  system. 
From  this  work,  we  can  conclude  that: 

(1)  The  PV-TE  hybrid  system  could  reach  the  highest 
efficiency  as  cutoff- wavelength  varies.  The  optimized 
cutoff  wavelength  of  the  system  is  mainly  determined 
by  band  gap  of  solar  cell. 

(2)  Cooling  systems  have  a  fundamental  influence  on  the 
performance  of  the  hybrid  system,  but  it  is  not  necessary 
to  use  heat  sinks  of  extremely  high  heat  transfer  coeffi¬ 
cients.  Proper  range  of  heat  transfer  coefficient  can  be 
figured  out  by  analyzing  the  system  performance. 

(3)  The  TEG  subsystem  plays  an  important  role  in  power 
generation;  it  contributes  about  10%  output  power  in 
a  GaAs-CoSb3  hybrid  system. 


(4)  Compared  with  PV-only  systems,  PV-TE  hybrid 
systems  are  more  suitable  for  working  under  high 
concentration  conditions.  Efficiency  of  the  solar  cells 
is  much  higher  in  PV-TE  hybrid  system. 

(5)  The  concentration  ratio  CJ  increases  almost  line- 
arly  with  heat  transfer  coefficient  of  the  cooling  sys¬ 
tem.  Different  values  of  y  indicate  the  feasibility  of 
the  system  working  under  high  concentrations. 

(6)  For  designing  the  PV-TE  hybrid  system,  to  increase 
concentration  ratio  and  use  a  suitable  heat  sink  simulta¬ 
neously  is  a  potential  way  to  improve  the  system  efficiency. 
Guidelines  for  system  design  have  been  provided. 
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Appendix  A. 

A.l.  Numerical  model  of  the  TE  generator 

The  TEG  subsystem  consists  of  a  heat  collector  with  a 
high-absorption  material  coated  surface,  a  thermoelectric 
(TE)  module  with  n-type  and  p-type  thermoelectric  mate¬ 
rial  legs,  and  a  heat  sink.  Fig.  13  shows  a  schematic 
diagram  of  the  solar  thermoelectric  generator,  n-type  and 
p-type  legs  of  the  TE  module  is  of  length  L  and  divided  into 
n  subsections  as  suggested  by  Li  et  al.  (2010).  The  advan¬ 
tage  of  this  finite  element  model  is  that  it  can  be  used  to 
estimate  the  performance  of  the  TE  module  accurately 
for  much  larger  temperature  gradients.  The  value  of  finite 
elements  n  should  be  large  enough  in  order  to  increase 
the  calculation  accuracy  and  obtain  a  constant  efficiency 
of  the  TE  module. 

For  subsection  layer  i  of  the  n-type  and  p-type  legs,  elec¬ 
trical  resistance  R^p,  thermal  conductivity  Kp,n  and  Seebeck 
coefficient  odN  can  be  obtained  as, 

I?*  = - IE -  (i4) 

'  A**  ■(,Gn’P{Tnifx)  +  Gn’P{Tni'p))/2  K  ’ 

Anp  •  +  Kn'P{T’’’p))/2 

L/n  [  ’ 

ot"’p  =  (<xnp(T’’fl)  +  «r’p{Tni’p))/2  (16) 

where  Tpf  is  the  interface  temperature  between  i  and  i  +  1 

subsections  of  the  TE  module,  An'p  is  the  cross-section  area 
of  the  n-type  and  p-type  legs;  onp,  Kn,p,  (xn'p  are  temperature- 
dependent  electrical  conductivity,  thermal  conductivity  and 
Seebeck  coefficient  of  n-type  and  p-type  legs. 

For  the  TE  module,  thermal  energy  transfers  between 
subsections  of  n-type  and  p-type  legs,  part  of  the  thermal 
energy  converted  to  electricity,  and  joule  heat  generated. 
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For  the  i- th  single  finite  element  of  n-type  leg  and  p-type 
leg,  the  absorbed  heat  is  (Wu,  1996), 


SS5  =  "(7T i  -  Tpn  -  -2i2teR7  +  *71teT7  (17) 


And  the  heat  remove  from  this  element  is 


07i  =  *r(7Ti  -  T7)  +)i2teR7  +  *71  teTp7 


(18) 


For  the  1st  element  of  the  TE  module,  the  heat  absorbed 
by  the  p-type  and  n-type  legs  equals  to  the  heat  absorbed 
by  the  hot-junction  of  the  TE  module, 


Qh,te  =  GTh,  1  +  GTh,  1 


(19) 


And  it  is  equal  to  the  heat  transferred  by  heat  collector, 

Qh.TE  =  Khc(Thc  -  Th)  (20) 


For  the  n-th  element  of  TE  module,  the  removed  heat  of 
p-type  leg  and  n-type  leg  equals  to  the  heat  removed  by  the 
cooling  system, 

Qcool,TE  =  Qc,n  +  0C,«  (21) 


I  f  Qconv^V 

|  PV  ^  QradPV 


Solar  Cell 


Qcool.PV 


T  coolant 


•|  Load  [*■ 


Fig.  14.  Schematic  representation  of  a  thermal  module  of  CPV. 


According  to  the  relationship  between  spectral  open- 
circuit  voltage  and  band  gap,  the  spectral  open-circuit 
voltage  of  a  solar  cell  under  1  x  SUN  can  be  expressed  as, 

>Wl)=^-&  (29) 

As  £jg 


And  the  heat  removed  from  the  i-th  element  is  the  heat 
absorbed  by  the  (i  +  l)-th  element  of  the  TE  module, 

=  &H%  (22) 


The  current  supplied  by  the  Seebeck  voltage  is  calcu¬ 
lated  as, 


Im  — 


V  te 

Rte  +  Rl 


(23) 


in  which,  the  expressions  of  Seebeck’s  voltage  Vje  and  total 
internal  electrical  resistance  RTE  are, 


n 


Y57U  -  T>"  +  7U  -  TpM] 

i=l 

(24) 

i=\ 

(25) 

where  Voc  is  the  open-circuit  voltage  of  the  solar  cell  under 
1  x  SUN  and  full-spectrum  conditions. 

Spectral  efficiency  of  solar  cell  under  1  x  SUN  is, 


h  pv,Xs(\) 


Foc^(1)-^(1)-FF 
Epv{ 1) 


(30) 


where  FF  is  the  fill  factor  of  the  solar  cell. 

As  the  concentration  level  is  increased,  the  spectral 
short-circuit  current  of  solar  cell  Isc,ki(Cg)  increases  linearly 
with  concentration  ratio  Cg  (Schwartz,  1982), 

ISc,is(Cg)  =  CgIsc^s{  1)  (31) 

and  the  spectral  open-circuit  voltage  Vocy{Cg)  increases 
logarithmically  with  Cg, 

VocACg)  =  vocM  1)  +  -  :  -  :  —  ■  In  Cg  (32) 
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The  maximum  output  power  of  the  TEG  subsystem  can 
be  achieved  when  RL  =  RTE, 

Rte  —  ( Ite)2Rl  =  {Ite)2Rte  (26) 


From  the  energy  equations  analyzed  before,  the  effi¬ 
ciency  of  the  TEG  subsystem  is, 


(27) 


A.2.  Numerical  model  of  the  solar  cell 


Fig.  14  shows  a  schematic  diagram  of  the  CPV  subsys¬ 
tem,  under  1  x  SUN,  the  spectral  short-circuit  current  den¬ 
sity  of  solar  cells  can  be  expressed  as  (Gopal  et  al.,  2005) 


Js 


C,is(l)  =  f 
J  21 


Is 

280 nm 


q  •  A  •  EQE  •  ■  F{X)  ^ 

hP  •  c 


(28) 


where  nf  is  the  diode  ideality  factor,  which  varies  from  1  to 
2  depending  on  the  manufacturing  process  and  semicon¬ 
ductor  material. 

Taking  the  influence  of  series  resistance  Rs ,  concentra¬ 
tion  ratio  Cg  and  cell  temperature  TPV  into  consideration, 
the  expression  of  the  CPV  subsystem  efficiency  must  be 
changed  to  the  following  equation  (Cheknane  et  al.,  2006), 

hpv  =  VlpvjA^)  +  Pn  ‘  {Rpv  ~  298)]  •  {1 

rif  •  kB  •  TPV  n 

+  q  *  [Foe, 2^(1)  +  PVoc  •  {Tpv  —  298)]  '  n  a 
•  [/^(l)  +  PIsc  •  ( TPV  -  298)]2  •  Rs 

77 

EpV 

where  ftIsc  and  PVqc  are  temperature  coefficients  of  the 
solar  cell. 

Output  power  of  the  solar  cell  can  be  written  as, 

Rpv  —  Opv  ‘  Rpv  (34) 
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Fig.  15.  Simulation  of  GaAs  solar  cell  and  comparison  with  experimental 
results,  Cg=  1,  TPV=  298  K. 

A  comparison  (Fig.  15)  of  the  solar  cell  efficiency  esti¬ 
mated  from  numerical  models  shows  good  agreement  with 
C.  Algora’s  GaAs  solar  cell  tests  (Algora  et  al.,  2001)  when 
the  concentration  ratio  is  ranged  from  1  to  2000  suns. 
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